radioactivedecay; A Python package for radioactive decay calculations

Malins, A.; Lemoine, T.*

Journal of Open Source Software (Internet), 7(71), p.3318_1 - 3318_6, 2022/03

https://doi.org/10.21105/joss.03318

radioactivedecay is a Python package for radioactive decay modelling. It contains functions to fetch decay data, define inventories of nuclides and perform decay calculations. The default nuclear decay dataset supplied with radioactivedecay is based on ICRP Publication 107, which covers 1252 radioisotopes of 97 elements. The code calculates an analytical solution to a matrix form of the decay chain differential equations using double or higher precision numerical operations. There are visualization functions for drawing decay chain diagrams and plotting activity decay curves.

Development of a code system DEURACS for theoretical analysis and prediction of deuteron-induced reactions

中山 梓介; 河野 広*; 渡辺 幸信*; 岩本 修; Ye, T.*; 緒方 一介*

EPJ Web of Conferences, 146, p.12025_1 - 12025_4, 2017/09

https://doi.org/10.1051/epjconf/201714612025

近年、重陽子加速器を使った大強度中性子源が様々な応用に対して提案されている。このような中性子源の設計のためには、高精度かつ広範な重陽子核データライブラリが必要不可欠である。それゆえ、我々はDEURACS (DEUteron-induced Reaction Analysis Code System)という、重陽子入射反応の解析と理論予測用の統合コードシステムを開発している。本研究では、反応の解析をこれまでよりも高い100MeV近傍まで拡張した上、DEURACSを80および100MeVにおける反応にも適用した。その結果、DEURACSの計算値はおよび反応の二重微分断面積の実験値を良く再現した。

国産安全解析コードの現状と課題; 我が国の安全規制への貢献を目指して

須山 賢也; 平尾 好弘*; 坂本 浩紀*

日本原子力学会誌ATOMO, 57(12), p.787 - 791, 2015/12

原子力施設等に対する世界最高水準の安全性を追求する取り組みにおいては、安全評価に使用する解析コード等の水準を高度に保ち続けることが求められる。我が国においては、それらの多くが原子力技術導入時に米国から導入されため、国産の解析コード等の開発が継続せず、結果的に古い米国産コードを利用しつづけている状況にある。本解説では、放射線工学分野における安全解析コードの現状と課題を示し、最新の知見を取り込んだ国産解析コードが我が国の安全規制や現場で広く利用されるための方策について述べる。

平成16年度大型計算機システム利用による研究成果報告集

情報システム管理課

JAERI-Review 2005-032, 151 Pages, 2005/08

JAERI-Review-2005-032.pdf:11.51MB

日本原子力研究所計算科学技術推進センターでは、研究活動を支援するためにスーパーコンピュータをはじめとする大型計算機システムを導入し、所内における研究活動に役立たせるとともに、計算機システム及びネットワークシステムの運用管理を行っている。本報告書は、平成16年度における原研の大型計算機の利用実績を集計し、利用ユーザからの利用報告に基いた研究内容,利用及びその成果についてまとめたものである。

モンテカルロ法臨界計算収束性ガイド資料; 原子力コード評価専門部会平成13,14年度活動成果

原子力コード研究委員会原子力コード評価専門部会

JAERI-Tech 2003-078, 107 Pages, 2003/11

JAERI-Tech-2003-078.pdf:6.18MB

核分裂性物質を含むユニット間の中性子結合の弱い相互干渉系を対象とするモンテカルロ法臨界計算では、解の収束緩慢性が問題になることがある。本報告では、この種の問題に対する適切な判断と解決方法の手がかりを与えられるように、基礎的な臨界計算の理論と収束性にかかわる応用理論を主として述べてある。このためには、この科学技術分野にかかわる論文の現状調査を広く行って、その内容を検討・参考とし、必要に応じてガイド資料に引用した。さらに、核分裂マトリックス固有ベクトルを利用したモンテカルロ法臨界計算の収束性加速及び判定方法の開発を行い、経済協力開発機構原子力科学委員会(OECD/NEA/NSC)の臨界収束性専門家会合で提案されたベンチマーク問題の多種多様な物理体系により手法の検証を行った。

弱結合相互干渉系臨界解析のマトリックス固有ベクトルを用いた収束性加速

野村 靖; 高田 友幸; 角谷 浩享*; 黒石 武

JAERI-Tech 2003-020, 88 Pages, 2003/03

JAERI-Tech-2003-020.pdf:4.31MB

弱結合中性子相互干渉系の中性子増倍率をモンテカルロ法計算コードにより求める場合に、しばしば解の収束性に関する問題が生じる。本報告では、この計算コードによる解の収束性に関して、一般的な中性子輸送方程式から導いた理論により、マトリクスK手法による解の収束性加速について考察を加えた。これにより、連続エネルギーモンテカルロ法コードMCNPにマトリクスK計算機能及び収束加速手順を組み込んだ。さらに、新たに開発した計算コードを用いてOECD/NEA臨界計算ベンチマークの2問題について解析を試み、マトリクスKによる解の収束性加速の有効性を確認した。

Marine environmental assessment system of radionuclides around Japan

小林 卓也; 外川 織彦

Proceedings from the International Conference on Radioactivity in the Environment (CD-ROM), 4 Pages, 2002/09

海洋中へ放射性物質が放出される際に短期間(30日)の核種移行予測と日本人の集団線量を算出する海洋環境影響評価システムSTEAMERを開発中である。短期間の核種移行予測コードシステムは海流場を解析するプリンストンオーシャンモデル(POM)と溶存放射性核種の海洋中拡散を解析するランダムウォークモデルSEA-GEARNから構成される。予測結果の可視化システムも同時に開発した。本システムの適用試験として、海流,水温,塩分,領域客観解析気象データ(RANAL)を用いて原子力潜水艦の仮想沈没事故の影響解析を実施した。核種拡散結果を比較するため、そして日本人に対する集団線量を算出するために、もう一つの計算コードDSOCEANも使用した。DSOCEANは等分割コンパートメントモデルと地衡流解析モデルから構成される。

Monte Carlo 2000国際会議の報告

中川 正幸

炉物理の研究, (51), p.47 - 51, 2000/12

http://rpg.jaea.go.jp/else/rpd/annual_report/UptoNo53/No51.pdf

MC2000(Advanced Monte Carlo on Radiation Physics, Particle Transort Simulation and Applications)国際会議が平成12年10月23日から26日までポルトガルの首都リスボンで開催された。筆者は本会議のgeneral chairmanとして会議の準備の段階からかかわっていたので、この会議の報告について執筆を依頼された。会議の参加者が303名,参加国30か国,総発表件数は237件であった。日本からは約20名が参加された。この会議はタイトルにあるようにモンテカルロ法を用いて放射線粒子の反応と輸送を扱う分野に限定していたので、すべての発表論文は何らかの形でモンテカルロ法と関係があるものばかりである。すべての放射線粒子が対象といっても、その大部分は(1)電子-光子系,(2)中性子-線,(3)ハドロン粒子系に分けられる。したがって会議もこの3分野を柱として構成されている。実際に発表された論文は(1)が114件,(2)が70件,(3)が50件であった。これに加えるに初日の全体会合で各分野から1件ずつの基調講演が行われた。この報告では、本紙が原子力学会の炉物理部会の会報であることから、中性子-線セッションを中心に幾つかの話題について紹介する。

「モンテカルロ計算夏季セミナー」実施報告

桜井 淳; 久米 悦雄; 谷田部 茂*; 前川 藤夫; 山本 俊弘; 長家 康展; 森 貴正; 植木 紘太郎*; 内藤 俶孝*

日本原子力学会誌, 42(10), p.1062 - 1065, 2000/10

日本原子力学会「モンテカルロ法による粒子シミュレーション」研究専門委員会主催の「モンテカルロ計算夏季セミナー」を実施した。大学・研究機関・企業から計111名の参加者があった。本稿においては、開催主旨、ノート型パソコンへの連結エネルギーモンテカルロ計算コードのインストールの技術的問題検討、セミナープログラム、計算メカ、参加者からのアンケート、今後の課題についてまとめた。今回は最初の読みであったため、実施上のさまざまな問題が生じたが、今後のために役立てていきたい。

ナトリウム漏洩燃焼形態の予測手法に関する研究 -粒子法を用いたナトリウム漏洩燃焼解析手法の開発- 中間報告書

越塚 誠一*; 向原 民*; 岡野 靖; 飯田 将雄*; 山口 彰

JNC TY9400 2000-012, 91 Pages, 2000/03

JNC-TY9400-2000-012.pdf:2.82MB

液体ナトリウムの漏洩燃焼挙動の解析には、ナトリウム液滴やナトリウムプールの燃焼、床に落下したナトリウムの広がり、燃焼生成物の堆積など、様々な現象を複合して解析する必要がある。特に、床ライナーの温度分布の評価には、落下したナトリウムがどのように広がるかが重要である。核燃料サイクル開発機構では、こうした複雑な液体ナトリウム漏洩燃焼挙動解析のため、粒子法による計算コードを開発している。粒子による熱流動解析手法は東京大学において本研究者らによって開発されたもので、従来の差分法などと比較して、流体の分裂・合体や堆積・凝固を容易に扱うことが可能である。本研究では、ナトリウムの凝固や燃焼生成物の堆積挙動を解析するため、固相の粒子計算モデルの開発、これを利用した凝固・堆積の粒子計算モデルの開発、表面張力の粒子計算モデルの開発をおこなった。固相の解析のため、ヤング率とポアッソン比で記述される弾性体の運動を、粒子間相互作用で計算するモデルを新たに開発した。この方法では、従来の有限要素法での要素分割の必要がないので、弾性体の大変形のみならず、破壊なども容易に解析できる。特に粒子に回転の自由度を持たせたことで、角運動量の保存性が得られた。弾性体に正弦波状の変位を与え本計算モデルを適用したところ、応力分布や圧力分布は解析解と良く一致した。凝固を伴う熱流動解析法として、熱流動をこれまでの粒子法(MPS法)で、凝固した後の固相を弾性体として解析するモデルを開発した。初期に幅10cm高さ20cmの2次元矩形状液体ナトリウムの、厚さ1cm長さ1mのステンレス板上での広がり挙動を、本手法を用いて計算した。仮にナトリウムおよびステンレスの熱伝導率や粘性を実際よりも大きな値を用いると、先端が凝固することにより、ナトリウムの広がりが抑制されるという結果が得られた。なお、本計算モデルは燃焼生成物の堆積挙動にも適用できる。表面張力の効果を取り入れるため、その粒子計算モデルを開発した。液面形状を描かずに粒子数密度から曲率など必要な微分幾何の諸量を計算するアルゴリズムとし、流体が分裂や合体をする場合にも適用できるようにした。エタノール液滴の振動の計算をおこない、従来のVOF法による結果と良い一致を得、モデルが妥当であることを示した。さらに、表面張力が卓越する場合の流体の広がり挙動を解析し、広がった流体が複数の液滴